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Abstract 



Novel systems of layered graphene are attracting interest for theories and applications. The stability, band structures of few-layer graphite 
films, and their dependence on electric field applied along the c-axis are examined within the density functional theory. We predict that those of 
Bernal type and also rhombohedral type tri- and tetra-layer graphite films exhibit stability. The rhombohedral-type systems including AB-bilayer, 
show variable band gap induced by perpendicular electric field, whereas the other systems such as the Bernal-type films stay semi-metallic. 
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1. Introduction 

Amongst the variety of allotropes of carbon over the full 
range of dimensionality, two-dimensional (2D) crystals consist 
of a single or few graphite layers have been stimulating sci- 
entist's interest in theories and in possibility for novel appli- 
cations. The observations of anomalous electric field effect on 
Hall coefficientlT] , unconventional integer quantum Hall effects 
in monolayerfSl and bilayer graphitefT] have been enabled by 
recent achievement of high quality preparation of the ultra-thin 
graphitic films IT]. 

These unusual electronic properties are believed to originate 
in unusual band structures near Fermi level of graphene[4J, a 
graphitic monolayer, and of its layered structures. Graphene 
has a direct zero energy gap between fc-linear tt bands crossing 
at each corner (K point) of the 2D hexagonal Brillouin zone. 
■^The charge carriers in graphene therefore are massless Dirac 
fermions, which leads to shifted Hall plateaus[2|. While the chi- 
ral parabolic band dispersion in AB-stacked bilayer graphene 
was found to give another new type of Hall effect||3|. 

Very recently several groups have reported calculated low- 
energy band structures of AB-stacked graphene layers using 
Tight-Binding (TB) models for tt band at zero electric field|5| 
and with field perpendicular to the layer|61. In the latter, they 
examined field induced band gap in AB-bilayer and found inter- 
esting band deformation in ABA- and ABAB-stacked systems. 

However, reliable TB parameters are only available for 
the limited equilibrium structures. Furthermore, the self- 
consistency should be crucial for studying their response 
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against the electric field, because of the induced charge re- 
distributions or polarizations as pointed out in the studies of 
dielectric function of graphite|7|. Properties of non-AB-type 
systems must be worth examining, because 16 percent of the 
natural and synthesized graphite are rhombohedral (ABC- 
stacked) with other 4 percent disordered. 

The purpose of the present study is to examine the band 
structures and interlayer binding energy of few-layer graphite 
(FLG) within the density functional theory (DFT) with the lo- 
cal density approximation (LDA), focusing on the effects of 
stacking sequence and external electric field. 

In the next section, we summarize the method of calculation, 
then possible equilibrium structures are sought in Sec. [3] Sec- 
tion |4] illustrates the calculated band structures near the Fermi 
level, followed by a concluding section. 



2. Methods 

We employed the CPMD pseudopotential planewave codelU 
as the ab initio total energy and band calculations in the present 
study. Goedecker-Teter-Hutter pseudopotential for carbon to- 
gether with Fade approximant LDA exchange-correlation po- 
tential were usedf9l. A sufficiently large planewave cutoff of 
120 Ry was used to achieve convergence for the wavefunctions, 
and 480 Ry for the density. The Monkhorst-Pack lfTOl special 
fc-point mesh of 10 x 10 x 1 within the first Brillouin zone (BZ) 
was used in the density optimizations. 

Two-dimensional (2D) structures of two up to four graphene 
layers are simulated using 3D hexagonal unit cells defined by 
primitive lattice vectors (see Fig. [TJ ai — a(l,0,0), a2 = 




Fig. 1. An in-plane translation (AtoB) of an upper graphene layer relative 
to the bottom layer. Triangle ABC defines the periphery of irreducible trans- 
lations. 

a(— i, ^\/3, 0), and as = (0, 0, c) with the optimized param- 
eter a ^ 2.4595 A. and c = 15, 20 and 25 A for bi-, tri- and 
tetra-layer systems, respectively, which ensure at least 10 A of 
vacuum separations along c-axis. 

The primitive reciprocal vectors in the basal plane are bi = 
(27r/a)(l,l/\/3), b2 = (27r/a)(0, 2/\/3), and the corners of 
the first Brillouin zone are labeled "K" which are given by 
i(bi + b2) plus the arbitrary reciprocal vectors. 

We examine the stability of FLG systems using the total en- 
ergy calculation within the DFT-LDA. It is known that DFT- 
LDA fails to describe long-range dispersion (van der Waals) 
interactions, but it well describes the short-range repulsion fTTj 
and predicts the ground state structure and interlayer binding 
energy [12 | in good agreement with experiments. Furthermore, 
Lennard- Jones type van der Waals energy does not make much 
difference between a discreat-atom and continuum models for 
graphitic plane. The zero-point phonon energy may not be im- 
portant in the present calculations, since the influence of stack- 
ing configuration on it is considered to be negligible. 

3. Structures of FLGs 

Figure|2a) shows the calculated total energies of AX-bilayer 
graphite with different stacking configurations compared to that 
of the Bernal AB-bilayer, in which the translation of the top 
layer X is B (and equivalently C). The interlayer separation 
(to be referred to as d) has been optimized for each structures. 
Note that energies are represented in meV per carbon atom on 
the top overlayer, which in fact correspond to the exfoliation 
energies 1, 12J . We see that AA structure marks, as is well-known, 
the highest energy of about 10 meV/atom with d = 3.62 A 
compared to the lowest energy of AB (AC) structure with d = 
3.33 A. This energy of AA-bilayer is as large as a half of the 
exfoliation energy observed (22.8 meV/atom) or predicted (24 
meV/atom) for graphite | IT2l . A translation of the top layer X 
from the in-plane position B to the nearest minimum C will 
encounter adiabatic energy barrier of 1.82 meV/atom at the 
midpoint, which is about 20 percent of the energy increase at 
the AA stacking. 
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Fig. 2. Energy surface of FLGs, (a) AX-bilayer, (b) ABX-trilayer, (c) 
ABAX-tetralayer and (d) ABCX-tetralayer, as a function of in-plane displace- 
ment of the top graphene overlayer (X). The total energies of the Bemal-type 
(AB, ABA and ABAB) stacking are taken as reference energies. Note that 
energies are represented in meV per carbon atom on the top layer. 

This kind of barrier height is more pronounced (2.92 
meV/atom) in ABX-trilayer graphite as shown in Fig. EJb), in 
which the structure of the first two graphene layers are fixed at 
the AB-bilayer optimized by itself and the top layer is allowed 
to translate. The separations from the substrate AB-bilayer are 
optimized. ABC-trilayer with optimum separation of d — 3.33 
A is found to have a slightly higher energy (0.18 meV/atom) 
than ABA with d = 3.33 A. 

The behavior of barrier along the translational paths 
(Figs.|2c) and (d)) is interesting in tetralayer graphite. The path 
linking ABAC (0.18 meV/atom) and ABAB (0 meV/atom by 
definition) has only a reduced barrier peak of 1 .27 meV/atom, 
which would make transformation to Bernal stacking easier, 
whereas the path from ABCA (0.33 meV/atom) to ABCB 
(0.18 meV/atom) has a much harder barrier peak of 3.07 
meV/atom. Since ABCB is equivalent to ABAC, we find that 
the most stable Bernal-type ABAB -tetralayer and the second 
most stable rhombohedral-type ABCA-tetralayer is separated 
by a translational barrier of about 3 meV/atom, if only a sur- 
face layer is allowed to move. These stable structures have a 
common optimum interlayer separation of d = 3.33 A. This 
result is consistent with the observation of rhombohedral struc- 
ture in natural and synthesized graphite fl^, and possibility of 
fabricating non-Bernal tri- and tetra-layer graphene films can 
be speculated. We have examined the effect of external elec- 
tronic field on the relaxed structures of these FLGs. However, 
no discernible change in the order of stability and interlayer 
separations was found. 

4. Band structures and field effect 

For FLGs that have been found globally or locally minimum 
in total energy (i.e., AB-bilayer, ABA- and ABC-trilayer, and 
ABAB-, ABCA-, ABAC-, ABCB-tetralayer systems), we have 
calculated the band structures of the systems in the zero field 
and also in an uniform external electric field of strength £ = 
77.94 mV/A, applied perpendicular to the graphene plane. This 
field was represented by an additional sawtooth potential, which 



is piecewise proportional to z with discontinuity at the mid- 
plane of the vacuum region in the cell. The potential difference 
applied between two adjacent graphene planes is = 0.2595 
Volts. 

4. 1 . AB-bilayer graphite 

The band structure near the corner of BZ (K point) of AB- 
bilayer graphite with no external field is semi-metallic with 
chiral parabolic dispersions as shown in Fig.[3ja). The effect of 
perpendicular field £ opens an energy gap as shown in Fig. [3b). 
The gap at the K point may indeed be understood in terms of the 
TB description (such as the SWMcC model| 13|) of tt electrons 
as a lifted degeneracy of non-bonding pz orbitals, as a conse- 
quence of the effective field that breaks the symmetry along 
the c-axis. The TB Hamiltonian of AB-bilayer graphite at the 
K point is identical to that for the system consists of one dimer 
and two equivalently isolated atoms as depicted in Fig. Oc). 
Indices 1 and 2+'^^ in Figs. Oa) and (b) label monomeric 
non-bonding and dimeric bonding(antibonding) states at the K 
point. The hopping parameter for the dimer, which can be de- 
termined by a half the difference in eigenvalues of 2* states at 
zero field, is found to be 71 — 0.363 eV in agreement with the 
well-known value of the SWMcC model. 

We point out that the gap at the K point, E^^{K) = 0.117 
eV, is only 45 percent of the unscreened potential difference 
AV (indicated by the length of a bar in Fig. Oc)) that the 
non-self-consistent TB model should predict. Our calculated 
result of the gap directly measures the effective field felt by the 
AB-bilayer film, and hence the dielectric response of bilayer 
graphite film against uniform and static perpendicular field, 
giving the dielectric constant of e(AB) = 2.2. This value of the 
internal effective field was also confirmed from the calculated 
lowest two eigenvalues of a band at the BZ center(r point). 

In addition, no discernible influence of the field was found 
on the hopping parameter 71, because the field induced gap 
E^^(K) was perfectly consistent with that predicted by 71 at 
zero field and effective potential difference e^^AV, i.e., 2['yf + 
(e-iAy)V4]^/2. Moreover, this fact was further confirmed in 
a separate calculation, in which the strength of the external field 
was doubled to widen the band gap at the K point to 0.221 eV. 
Meanwhile, Ohta et al.flA] have reported significant change 
in the dimeric interaction with electron doping by potassium 
adsorption on a surface of AB-bilayer, from their angle-resolved 
photoemission experiments with analysis using TB fitting. This 
apparent contradiction with the present results might be related 
to the shift of Fermi level (electron number) due to doping in 
their system. 

4.2. Trilayer graphite 

Figure |4|a) shows that the band structure near the corner 
of BZ (K point) of ABA-trilayer graphite film at zero field 
is semi-metallic with rather small overlap of conduction and 
valence tt bands. The TB Hamiltonian of the ABA-trilayer at 
the K point is equivalent to that for a linear trimer(see panel 
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Fig. 3. Band structures of AB-bilayer grapliite film: (a) without field, (b) 
with field of £■ = 77.94 mV/A, (c) shows TB diagrams at the K point. The 
inset magnifies 1/8 range of (a) about K and energy range of ±30 meV. For 
indices 1 and 2^ see text. 
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Fig. 4. Band structures of trilayer graphite films: (a) ABA-stacked without 
field, (b) with field of £ = 77.94 mV/A, (d) ABC-stacked without, and (e) 
with the field, (c) and (f) are TB diagrams at the K point. The inset in (a) 
and (d) magnifies 1/8 range about the K and energy range of ±30 meV. 

(c)) with hopping for the nearest neighbors 71 and for the next 
neighbor 75 that links both ends, and a perfectly non-bonding 
monomer, and a very weakly bonded (with 72; a dashed link 
in panel (c)) dimer that opens a tiny gap (12 meV) between 
nearly fc-linear branches. This trimer at zero field has three 
states, which are to be referred to as bonding/antibonding and 
non-bonding states labeled by 3+^^ and 3", with eigenvalues 
(^75 =F ^^27^ + 7I /4) and —75, respectively. Note that 3° 
actually is the bonding state of orbitals at both ends via the 
second neighbor hopping 75 . 

On the other hand, ABC-trilayer graphite film at zero field is 
as shown in Fig. |4|d) a zero-gap semiconductor(or zero band- 
overlap semi-metal) with a band contact near the K point on the 
KM symmetry line. The gap at the K point is 20 meV, which 
arise from the splitting between bonding(l+) and antibonding 
(1~) states the weakly bonded dimer with the hopping param- 
eter of 72 = 10 meV. The strongly bonded dimers (Fig. 15 f)) 
are giving degenerate states. 

The influence of the electric field on band structures of tri- 
layer systems are shown in Fig.Ud) and (e). We see the ABA- 
trilayer stays semi-metalHc, because of the presence of almost 
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Fig. 5. Band structures of tetralayer graphite films, (a) ABAB-stacked without 
field, (b) with field of £ = 77.94 mV/A, (d) ABCA-stacked without, and 
(e) with the field, (g) ABAC-stacked without, and (h) with the field, (j) 
ABCB-stacked without, and (k) with the field, (c), (f) and (1) are TB diagrams 
at K point. 

degenerate states of non-bonding monomer in the middle plane 
(B) and trimer, which are nearly independent of the applied 
field. Whereas, the ABC-trilayer opens a energy gap at the K 
point, E^^^{K) = 0.188 eV, and a minimum direct gap 0. 146 
eV. This energy gap is the 1+-1^ splitting and is nearly dou- 
ble the width of AB-bilayer's induced by the same external 
field, reflecting the doubled thickness of the film, since the 72 
is order of magnitude smaller than the external potential differ- 
ence. Therefore, the ABC-trilayer can be a tunable narrow-gap 
semiconductor if fabricated. The dielectric constant of ABA 
and ABC trilayers at the present applied field may be evaluated 
from the 1+-1^ splitting, which results in e(ABA) = 2.4 and 
e(ABC) = 2.7. 

4.3. Tetralayer graphite 

The band structures near the K point of four types of stable 
and metastable tetralayer graphite films are shown in Fig. |5] 



ABAC- and ACBC-stacked systems are obviously equivalent 
at zero electric field and even at an uniform field along c- 
axis they are equivalent when the reversed field is applied in 
either system. Therefore, the panel (k) may be thought of as 
the counterpart of panel (h) with reversed field. 

The Bernal-type, ABAB -tetralayer system at the K point is 
equivalent to a linear tetramer and two weakly bonded dimers 
configured as is shown in panel (c). This clearly explains the 
distribution of Kohn-Sham eigenvalues at the K point in panel 
(a). The labels 4* refer to the full bonding and antibonding 
states, and 4** are semi-bonding and antibonding states of the 
linear tetramer The external field lifts the degeneracy of upper 
and lower dimers. However, this does not leads to a band gap 
as show in panel (b) at this level of the field strength. From the 
splitting of weakly bonded dimeric states, the dielectric constant 
in the system is evaluated to be e(ABAB) — 2.4, almost no 
change from the Bernal trilayer's value. 

Similarly to the trilayer case, the rhombohedral-type ABCA- 
tetralayer system develops a flat or rather concave band gap 
as an uniform field comes in as shown in panel (e). The en- 
ergy splitting of nearly non-bonding states (labeled as 1) on 
the top and bottom layers is 0.263 eV, which corresponds to 
e(ABCA) = 3.0. The induced narrow band gap will be con- 
trolled using a variable field. 

The ABAC- or ABCB-stacked system has a very narrow 
band gap structure as shown in panel (g), which is rather similar 
to the semi-metallic band of ABA-trilayer around the K point 
at zero field, except for an overlaying dimeric band. When an 
external field is applied parallel to or anti-parallel to c-axis, the 
system becomes semi-metallic as plotted in panels (h) and (k). 
The 3'-^ and 1+ states are always very close, nearly independent 
of the uniform field, since they feel nearly the same averaged 
external potential. The dielectric constant is evaluated to be 
e(ABAC) = 2.6 and e(ABCB) = 2.5 at the present finite 
field, where we notice very small asymmetry due to the broken 
symmetry of the stacking sequence along the c-axis. 

5. Conclusions 

We have examined the stability of FLGs and found that for 
any number of layers Bernal FLGs are the most stable, but 
rhombohedral-type tri- and tetra-layer graphite films are also 
predicted to be possible structures, which are stabilized by rel- 
atively high energy barrier. In the rhombohedral-type FLGs in- 
cluding AB-bilayer, variable band gap will be induced by the 
on-set of perpendicular electric field, and other FLGs will stay 
semi-metallic. (Quite recentry, a report on ab initio band struc- 
tures of FLGs at zero fieldfTS'l came into our notice. We found 
that the present results at zero field is in fine agreement with 
their results.) Therefore, experimental challenges for the fabri- 
cation of these non-Bernal FLGs would be very interesting. 
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